Persistent current in metals with a large dephasing rate 
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In a weakly disordered metal electron interactions are responsible for both decoherence of the 
quasi-particles as well as for quantum corrections to thermodynamic properties. We consider elec- 
trons which are interacting with two-level-systems. We show that the two-level-systems enhance the 
average equilibrium ("persistent") current in an ensemble of mesoscopic rings. The result supports 
the recent suggestion that two puzzles in mesoscopic physics may be related: The low temperature 
saturation of the dephasing time and the high persistent current in rings. 
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Quantum interference effects play a crucial role in the 
low temperature properties of normal metals. Prominent 
examples are weak localization and the associated mag- 
netoresistance. Recently it was suggested [|]J^ that two 
of the unresolved problems in the physics of mesoscopic 
metals may have a common solution: The large value of 
the persistent current in mesoscopic rings and the low 
temperature saturation of the dephasing rate which is 
seen in magnetoresistance measurements. 

The first problem is the large value of the persistent 
current in rings. Levy et al |^ measured the nonlinear 
response to a magnetic field of an ensemble of 10^ meso- 
scopic copper rings. The measured signal corresponds 
to a current I k, Iq svn.{2iT(j)/ (jj^) circulating in each ring. 
(j) is the magnetic flux which penetrates each ring and 
(po = h/e is the flux quantum. For temperatures in the 
mK regime the amplitude was |/o| ~ 0.3nA per ring, 
which is of the order of one elementary charge in the 
time td an electron needs to diffuse around the ring, 
|/o| « 0.6e/TD = O.GeEjh. Here = H/td = hD/L^ is 
the Thouless energy, D is the electron diffusion constant, 
and L is the circumference of the ring. Similar results 
were reported in Refs. 0,^. 

Theory, when neglecting electron interactions, predicts 
a current that is of the order / ~ e6/h, where d is the av- 
erage distance of single particle levels at the Fermi energy 
With the parameters of the experiment ||], 6/k « 
0.2mK and Ec/k « 25mK, the current obtained is about 
two orders of magnitude too small. Electron interactions 
first seemed to improve the situation [^. For Coulomb 
interaction it was found that / ~ efi*Ec/h, where /i* is 
a dimensionless number that characterizes the strength 
of the interaction in the Cooper channel. However es- 
timates of /i* gave a value which is an order of magni- 
tude too small when comparing it with the experiment 

. Surprisingly an enhancement of the current was also 
reported in presence of a moderate concentration of mag- 
netic impurities with / ~ e{Ec/h) •min(?i/Ts, Ec)/kT, 
where l/r^ is the spin-flip scattering rate. This can be- 
come larger than the current coming from the Coulomb 
interaction, however since the temperature dependence is 



different from the one observed this mechanism has not 
been considered as a possible explanation of the experi- 
ment in Ref. 

The second problem concerns the phase coherence of 
the electrons. Whereas it is expected that the dephas- 
ing rate goes to zero in the zero temperature limit pC| ] 
many experiments show a saturation at low temperature. 
Usually this saturation is attributed to the presence of 
magnetic impurities or to heating. However, recently 
a saturation of the dephasing time has been observed, 
also after excluding these possibilities |l^,|T^. Several 
attempts have been made to explain the low tempera- 
ture saturation of the dephasing time [p"3|-pT|. It has 
been argued by Altshuler et al that non-equilibrium 
electromagnetic noise can decohere the electrons without 
heating them. Originally, this non-equilibrium noise was 
suggested to be due to external radiation which couples 
into the samples. On the other hand dephasing could 
also occur due to internal noise. In this case a saturation 
of the dephasing time could also occur in equilibrium. 
Experimental evidence is in favor of an internal dephas- 
ing mechanism ||ll|,|l^ , however it is open if equilibrium 
or non-equilibrium processes dominate. 

Recently Kravtsov and Altshuler ^ have extended ear- 
lier work on the effect of a high frequency electro- 
magnetic field in mesoscopic rings and have shown that 
non-equilibrium noise leads to a directed non-equilibrium 
current. They then suggested that both the "large" cur- 
rents observed in 0J^,01 and the strong dephasing are 
related and non-equilibrium phenomena. 

In this paper we demonstrate that also for the system 
in thermal equilibrium an enhanced persistent current 
is expected if there is an additional electron interaction 
which gives also rise to strong dephasing. For the partic- 
ular model involving two-lcvel-systems (TLS) we find (1) 
a diamagnetic current in the low magnetic field limit (2) 
a temperature dependence which is close to the experi- 
mentally observed one (3) an amplitude which depends 
on the concentration of TLS. In the following we first 
recall some of the theoretical concepts concerning the 
persistent currents. We then estimate the persistent cur- 
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rent coming from TLS and, finally, relate the persistent 
current amplitude and the dephasing rate. 

The equilibrium current in a ring which is penetrated 
by a magnetic flux is calculated by taking the deriva- 
tive of the thermodynamic potential, = — ^ri(/i, (f>). 
For simplicity we do not discuss the subtle questions con- 
cerning differences between the canonical F{N, (j)) and 
the grand canonical thermodynamical potential r2(/x, (j)) 
In an ensemble of weakly disordered rings the 
disorder configuration will change from ring to ring, so 
in order to calculate the average persistent current of 
an ensemble of rings one has to average over disorder, 
~^ (-^(0))dis- For non-interacting electrons the grand 
canonical potential, and therefore the persistent current, 
depends only exponentially weak on the magnetic flux 
and one finds only a small persistent current | p9[ . 

The situation changes in presence of interactions. As 
an example take the Coulomb interaction as in Ref. ||] 
and consider the classical expression for the Coulomb en- 
ergy. 



= - y drdr'w(r — r')Sn{r, (j))Sn{r' , cf)). 
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This quantity depends on the magnetic flux (p even on 
average since the fluctuations of the electron density are 
magnetic flux dependent and may be written as 

((5n(r, (^)fc(r', (^))dis = ^ cos(47rm0/(/)o) (2) 



with 

_ AN{eF)s\n\kp\v-v'\) ,^^^ 



(3) 



compare Ref. |8|j9|]. Here uj = 2nnkT are bosonic 
frequencies and V is the volume and N{eF) is the 
density of states of the Fermi level. As a tech- 
nical remark we would mention that Eq.(^ is ob- 
tained by evaluating the diagram with one particle- 
particle propagator (cooperon). The harmonics of 
the persistent current I{(f>) = /m sin(47r(/)TO/(/)o) 

are finally found as /,„ = IGfioe/Trm^TD with /iq = 
N {ep) J drv{r) sin^{kpr)/ (kpr)"^ . Including the ex- 
change energy reduces the current by a factor two, and 
higher orders in the interaction reduce the interaction 
amplitude, — > ~ A'o/il + Mo hi(eF'''D)]- 

When opening an additional interaction channel one 
will find an additional contribution to the persistent cur- 
rent. In Ref. this has been demonstrated for magnetic 
impurities. Here we consider the interaction of conduc- 
tion electrons with nonmagnetic impurities, which we as- 
sume to couple to the electron density. The Hamiltonian 
is of the form 



Hint = / dxn(x)l/(x). 



(4) 



The operator ^(x) that is due to the impurities will 
be specified more explicitly below. To second order in 
this interaction one finds a correction to the free en- 
ergy which is the sum of a Hartree and a Fock like term, 
6il ~ jr^H + Sflp , which are given by {P — 1/ kT) 
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At I dy. I dx'(n(x))(n(x')) 



(5) 



X [(f(x,r)f(x',0))-(f(x))(f(x')) 

X (vl/,(x, t)*^, (x', 0)) (f(x, t)V(x', 0)), (6) 

where 4'|(x, r) and ^'^(x, r) are operators for fermions 
with spin s and the brackets (. . .) are the thermal average 
. If F(x) describes pure potential scattering, then F(x) 
is a c-number with the result that JJIh = 0. 7^ but 
does not depend on the magnetic flux which can be traced 
back to the fact that ^(x, r) = V{'x) is static. This can 
become different if the impurity has an internal degree of 
freedom. Consider a TLS, realized by an impurity which 
sits in a double well potential with minima at r and r -I- d 
which are nearly degenerate in energy. We write the scat- 
tering potential as ^(x) = y[nA^(x— r)-|-ns(5(x— r — d)]. 
fiA and fiB are the number operators for the impurity in 
the relevant potential minimum. Since the impurity is in 
either of these minima fiA +nB = 1. We further charac- 
terize the impurity by the asymmetry e and a tunneling 
amplitude A between between the two minima, so the 
impurity Hamiltonian is 



e A 

A -e 



(7) 



The Hartree energy (||), which is nonzero in this model, 
may be interpreted from the point of view of both the 
electrons and the impurities. From the electronic point of 
view the electron impurity interaction gives rise to an ef- 
fective electron-electron interaction: Comparing Eqs.(|^) 
and (|) one realizes that the Coulomb interaction is re- 
placed by an effective interaction 



w(x — x') 



dr 



{(V^(x,r)f(x',0))-(V^(x))(y(x'))} (8) 



due to the defects. From the impurity point of view the 
coupling to the conduction electrons changes the level 
asymmetry. 



e A 

A -e 



+ V{n{r)) 
A 



A 

-e + V{n{r + d)) 



, (9) 



which then changes the free energy as given in Eq.(|^) to 
second order in V. 
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We discuss the persistent current first in the most sim- 
ple situation, where we neglect the tunnel splitting A. In 
this case V^(x, r) is static so (Jfip remains independent of 
magnetic flux, as in the case of "normal" disorder. Here 
and below we will therefore concentrate on the Hartree 
energy. Using the relation fiA + fiB = 1 and averaging 
over "normal" disorder we can rewrite the Hartree energy 
as 



{Snn)dis = -\V\^Sn^r,(f>))dis[l- 



sin^ (kpd) 



X / AT{{hA{T)nAm-{hA){nA)}. (10) 

If the TLS asymmetry is large, |e| > kT, then 
(nA(T)nA(O)) — (nA)(?^A) = and therefore (^ilH)dis — 0. 
For a TLS with a small asymmetry, |e| < kT one finds 
{fiA{T)hA{Q))-{nA){nA) = 1/4 so that (5f2H)dis 7^ and 
a persistent current results. From the integration over r 
it follows that the current coming from a single defect is 
proportional to the inverse temperature, in full analogy 
to the persistent current from a magnetic impurity [^. 
For the system with a finite density of TLS the asym- 
metry e will not be a constant, instead there will be a 
distribution of asymmetries. Using eq.(^) and below we 
determine the persistent current as 



kT td 



(11) 



where F = 1 — siv? {kpd)/{kpd)'^ and Cact is the concen- 
tration of TLS with e < kT and therefore is active in 
producing a persistent current. Assuming a fiat distri- 
bution of asymmetries between zero and Cmax > kT ^ the 
concentration of active TLS is proportional to the tem- 
perature, Cact = c/cT/emax, which then cancels the in- 
verse temperature dependence of the persistent current 
of a single defect. The current is diamagnetic in contrast 
to the paramagnetic current from the repulsive Coulomb 
interaction. The amplitude of the current is of the dif- 
fusive scale, / e/r/3, as for the Coulomb interaction. 
The dimensionless prefactor /i* is replaced by the factor 
A'TLS = —cFN{ep)V'^/e-a^ax- which should be of order 
one if this mechanism is relevant for the currents ob- 
served in Ref. Assuming an atomic scattering cross 
section of the TLS and the factor F ^ 1 this requires a 
density of states of TLS that is comparable to the density 
of states of the metallic host and therefore of the order 
10^^/Kcm'^. At lOOmK this corresponds to a concen- 
tration of active two-level-systems of about 2ppm which 
is not a small number but, in principle not impossible 
pO{ . For the assumed distribution of asymmetries the 
temperature dependence of the persistent current is only 
due to the temperature dependence of the local density 
fiuctuations, see eq.(|^), and is therefore identical to the 
temperature dependence of the persistent current from 
the Coulomb interaction. The latter has been shown 01 



to agree well with the experiment in Ref. Finally it 
is important to discuss spin-orbit scattering, since in the 
gold or copper rings in the experiments the spin-orbit 
rate is large. Following Refs. |8[|^ we find that spin- 
orbit scattering reduces the persistent current due to the 
mechanism discussed here by a factor four, but the sign 
remains diamagnetic. 

Let us now allow a finite tunnel splitting A, i.e. spon- 
taneous transitions of the impurity between the two min- 
ima. The correlation function that is relevant for the per- 
sistent current, i.e. the impurity susceptibility, is given 
by 



dT(ri^(T)n^(0)) - {nA){fiA) = 




4 e2+A2 Ve^+A^ 



(12) 



in the two hmits where -I- A^ < (fcT)^ and -I- A^ > 
{kT)^. Whereas for static defects with A = the corre- 
lation function is non-zero only in the high temperature 
limit, kT > e, the correlation function for dynamic de- 
fects is non-zero even in the zero temperature limit, so 
these defects contribute to the persistent current even 
for T 0. We calculate the persistent current under 
the assumption ||2l[ of a fiat distribution of e between 
zero and emax and a distribution of A that is propor- 
tional to 1/A between Amin and Amax- We then find 
/ ~ — (e/rD)F?i/(TTLsemax) as before when we neglected 
the tunnel splitting. ?i/ttls ~ cN{ep)V'^ is the electron 
scattering rate off the TLS. 

Finally we discuss the relation of the persistent current 
and dephasing. In Ref. |l^ it has been demonstrated that 
TLS lead to dephasing with a rate that is temperature 
independent in a certain range of temperature. Both the 
persistent current amplitude and the dephasing rate are 
hard to estimate for a given material since they depend 
on the concentration of TLS and the distribution of e and 
A. It is therefore of interest to relate the two quantities, 
in order to reduce the number of unknown parameters. 
Notice that in order to have dephasing there must be real 
transitions between two impurity states, and one finds 
that the defects with kT > yeM-A^ > Ti/r^ are most ef- 
fective for dephasing. On the other hand all defects with 
kT > vtM-A^ and even some with kT < Ve^ + A^ con- 
tribute to the persistent current, see Eq. (|T^) . We cannot 
therefore give a general relation between dephasing rate 
and persistent current amplitude. We can, however, as 
shown below, give such a relation for our special choice 
of the distribution of e and A. The dephasing rate has 
been estimated as [ITsI 
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max 



F/{e 



. max'^TLsA) 

max 

{F/hX 

^max 



LsA) if h/r^ < Amax < kT 

^max 

< h/T^ < kT 
(13) 
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with A = ln(Aniax/Amin)- The persistent current ampli- 
tude, / ~ MTLs(e/rD) with |/xtls| ~ Fh/ {erm,xTThs), is 
therefore of the order 



(14) 



in the two Hmits considered. For example for the gold 
sample of Ref. ^ fi/r^ ~ 2mK below 500mK. If the 
constant dephasing rate is from the mechanism we con- 
sider, then the lowest measured temperature (~ 40mK) 
is an upper limit for Amax, and leads to the estimate 
I/ZTLSI > A/20. 

The dephasing rate for low temperature, Ti/t^ < kT < 
Ainax, is proportional to T and given by l/r^ ^ 
Ffcr/(e,naxTTLsA). Here one finds |/iTLs| ~ X{h,/T^)/kT, 
which depends only on one unknown parameter, A. A 
dephasing rate which is linear in T has been observed in 
various three-dimensional and two-dimensional samples 
|p2| . The values which were reported for at lOK are 
around ^ lO^^^sec - 5-10~^°sec, which corresponds to 
{h/T(j,)/kT ~ 2 • 10^^ - 1. Also from these considerations 
it seems rather reasonable that the parameter |/iTLs| can 
reach values of order one. 

In this paper we estimate the persistent current lin- 
ear in the concentration of TLS and we neglect Kondo 
correlations. Kondo physics has been suggested as a pos- 
sible solution of the dephasing problem in Ref. |jl^ . The 
persistent current, of course, will be modified by Kondo 
correlations, however it is beyond the scope of this paper 
to estimate this quantitatively. We also do not attempt 
to give an exhaustive discussion of the limit of high con- 
centration of impurities. However it is clear that our 
theory will overestimate the current when the phase co- 
herence time becomes of the order of, or shorter than, 
the diffusive time tjj . The related problem for the persis- 
tent current coming from magnetic impurities has been 
discussed in Ref. 

In summary we demonstrated that the interaction of 
conduction electrons with impurities induces a persistent 
current. Under reasonable assumptions we find a tem- 
perature dependence that is set by the diffusive scale. 
The most crucial point however is the current amplitude 
here given by / ~ — |/iTLs|e/T£). The dimensionless pa- 
rameter ^TLS depends on the concentration of the TLS, 
so a reliable estimate of the current amplitude is dif- 
ficult. Experimentally the interrelationship of dephas- 
ing and persistent current may be checked by measuring 
the persistent current for different materials. For silver, 
where no saturation of the dephasing time has been ob- 
served [ p^ , we expect a smaller persistent current than 
in gold or copper where the dephasing time saturates at 
low temperature. The sign of the current may help to de- 
cide if non-equilibrium fluctuation suggested in Ref. |^] 
or the equilibrium electron-impurity interactions stud- 
ied here dominate the current: For a system with strong 



spin-orbit interactions Ref. |^ predicts a paramagnetic 
current, whereas we found a diamagnetic current. 
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